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We present results for the first two moments of the distribution amplitudes of pseudoscalar 
mesons. Using two flavors of non-perturbatively improved clover fermions and non-perturbative 
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compare with recent results from models and experiments. 
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1. Introduction 

Hadronic wave functions are of crucial importance when describing exclusive and semi-ex- 
clusive reactions O. For detailed references and applications consult pp. Distribution amplitudes 
(DAs) are related to the hadron's Bethe-Salpeter wave function, <j>(x,k±), by an integral over trans- 
verse momenta. For the leading twist meson-DAs we have 

<i>(x,p 2 )=Z 2 (p 2 ) / d 2 k(j>(x,k x ), (1.1) 

where x is the quark longitudinal momentum fraction, Z2 the renormalization factor (in the light- 
cone gauge) for the quark-field operators in the wave-function, and p denotes the renormalization 
scale. In this presentation we quote all numbers with a scale p 2 = 4 GeV 2 in the MS-scheme. 

It is convenient to rescale £ = 2x — 1. In the following we use 0(£) to describe any pseu- 
doscalar meson, to refer to the pion, and 0jr(£) to denote the kaon. Furthermore, it is 

common to expand DAs into their Gegenbauer moments and quote the expansion coefficients, a,-, 
at a given renormalization scale as a parameterization of DAs, 




(1.2) 



The zeroth moment is normalized to unity, f^d^Q^ ,p, 2 ) = 1, at any energy scale p 2 . From 
renormalization group arguments we find that 

Taking the u- and <i-quarks to be degenerate, G-parity implies that the pion DA is an even function 
of § and hence all odd moments vanish, i.e., afn+i = 0. 

Recently, we have computed the first moments of meson distribution amplitudes [Q] on the 
lattice. Independently, a calculation of the first moment of the kaon distribution amplitude has 
appeared which uses a different discretization scheme and different working points In this 
contribution we present our calculation of the first non- vanishing moment of the pion DA, aj, and 
the first two moments of the kaon DA, af and a\ . We compare our results to previous estimates 
from sum rules and experiment and discuss the implications of our lattice computation. 



2. Lattice calculation 

On the lattice, one has to perform the light-cone operator product expansion to find a relation 
between local operators and moments of DAs w.r.t. To be specific, one has 

<n(M 2 ) = J^dtrmn 2 ), 

(n|^ {vo ... v „ } (0)|PS) = if PsP{V0 ...p Vii} (^), (2.1) 
*V.v„ (0) = i" q(0)y v 75 d Vi ■ ■ ■ B v „ B(0) , 

where "PS" refers to either the pion or the kaon, and q can be either a d- or an s-quark D is the 
covariant derivative and {. . . } denotes the symmetrization of indices and the subtraction of traces. 
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The matrix element, Eq. (2.1), can be obtained from an appropriate ratio of two-point func- 
tions, consult [§] for details. For the first moment of pseudoscalar mesons containing non-degene- 
rate mass quarks we use the following two operators 

6%\ = ^{41} , < = ^{44} " ~ (^{11} + ^{22} + ^{33}) • (2.2) 

For we have used the external momentum p = (2jt/L,0,0) — with L being the spatial lattice 
size — and the corresponding rotated momenta, for we have taken p = 0. For the second 
moment we have only employed ^412 = ^{412} with external momentum p = (2n/L,27i/L,0). 

We have generated our gauge field configurations using the Wilson gauge action and two 
flavors of dynamical, non-perturbatively improved clover fermions. For four different values of 
j3 = 5.20, 5.25, 5.29, and 5.40 and up to four different K values per j8 we have produced ^(2000 — 
8000) trajectories. Lattice spacings and spatial volumes vary between 0.075 — 0.123 fm and (1.5 — 
2.2fm) 3 , respectively. The scale has been set using a Sommer parameter of ro = 0.467 fm. For 
further details see [Q]. Correlation functions are calculated every 10 HMC trajectories using four 
different locations of the fermion source. By applying binning we obtain an effective distance of 
20 trajectories. We observe that the bin size has little effect on the error, indicating that residual 
autocorrelations are small. 

The matching between the lattice results and the MS-scheme has been done non-perturbatively. 
The renormalization procedure has been detailed in Refs. [|5| |6|] and will be discussed further in a 
forthcoming publication. 



3. Numerical results 



3.1 Mass-degenerate quarks 

For mass-degenerate quarks, i.e., the pion, we compute (^ 2 )^(/x 2 = 4GeV 2 ), and from this 
«2 (4GeV 2 ). To obtain the result at the physical pion mass we first perform a linear chiral extrapo- 
lation of (^ 2 ) to the physical pion mass at each fixed value of j8. Figure [I] shows (i§ 2 ) obtained from 
^412 as a function of the pion mass at a fixed value of j8 = 5.29 with K" sca = K^. The dependence on 
the pion mass turns out to be very weak. Using chiral perturbation theory in the continuum it has 
been shown [{7|] that for small pion masses the leading logarithmic contribution can be absorbed 
in the pseudoscalar decay constant, f % . Hence, for small values of the pion mass we expect the 
dependence on m n to be rather flat. However, a precise matching of lattice results and chiral per- 
turbation theory similar to what has been done in [||] for the nucleon axial coupling, g&, is still to 
be performed. 

Figure || shows the continuum extrapolation of the results from all values of j3 at the physical 
pion mass for ^412- Scaling violations seem to be small, introducing an uncertainty of roughly 6%. 
The continuum result reads 

(^ 2 ) ;r ( A t 2 = 4GeV 2 )= 0.269(39), a%(4GeY 2 ) =0.201(114) . (3.1) 

It is in agreement with the older quenched result (<^ 2 )^(4GeV 2 ) = 0.286(49)1^ from [g]. It 
is larger than the asymptotic value, (^ 2 )^(/x 2 — > °°) = 0.2, indicating that the commonly adopted 
asymptotic ansatz at this energy scale may not be justified quantitatively. 
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Figure 1: Pion mass dependence of the second moment of the pion distribution amplitude, (^} 7r , at j3 
5.29, for all four values of fc sca = fc va] from ^412- 



The coefficient has been under experimental scrutiny recently. The decay of neutral pions 
into two photons contains information on the pion wavefunction. However, as has been reported in 
[10] the leading twist expression is insufficient to determine a? unambiguously. QCD sum rules, 



e.g. [11, O, 13] have found that 012 is indeed positive and the more recent calculations in this 



approach [ ]12| , |13J are compatible with our result also in magnitude. However, the contribution of 
aj is expected to be non-negligible in modeling the pion DA. This makes an independent lattice 
analysis of this quantity particularly interesting. 



3.2 Mass non-degenerate quarks 

In case of two quarks with distinct masses, the odd moments will no longer vanish. Since we 
are primarily interested in the light pseudoscalar mesons we tune the parameters and extrapolations 
in such a manner that the light quark (corresponding to fc seil ) has the correct mass to reproduce a 
pion if the quarks were degenerate and the heavy quark (corresponding to fe v m) then reproduces the 
kaon mass. Since we require a large set of different valence masses, we have restricted ourselves 
to the four data sets at j3 = 5.29 only. These have fc sea = 0.1340,0.1350,0.1355, and 0.1359. The 
systematic error due to the continuum extrapolation is estimated from the extrapolation in Fig. |^ to 
be about 6%. 

A sample plot for the second moment of the kaon distribution amplitude can be seen in Fig. || 
for the working point /3 = 5.29, fc; ea = 0.1359. Several different values have been chosen for the 
heavy quark mass. The vertical line indicates the physical kaon mass and the data point at this 
location has been taken from a linear extrapolation. 

The data from all four samples of gauge fields at j8 = 5.29 has been combined and linearly 
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Figure 2: Continuum extrapolation of (t; 2 ) K for all values of j3 for ^412. 

extrapolated to the physical values of and m n in a global fit, see for details. The resulting 
value is 



(% 2 ) K (4GeV 2 ) = 0.260(6)(16) , of (4GeV 2 ) = 0.175(18)(47) . 



(3.2) 



The first error is of statistical origin while the second one corresponds to a scaling error of 6%. The 
result (3.2) corresponds to a ratio of (^ 2 )k / (£, 2 )k — 1- I n the literature, there has been the debate 
between [11] predicting a ratio of af /a% ~ 0.59 ± 0.04 and [12, [H]] predicting a ratio of about ~ 1 . 
Our finding clearly favors the latter results. 

Of particular interest is the first moment of the kaon DA, and, hence, also af . Figure |] 
shows the interpolation of (E, ) obtained from as a function of the mass difference m\ — m 2 K for 
one sample of gauge fields. Note that this is one of the few cases where an interpolation suffices. We 
now have a result at a specific pion mass determined by k; m = JC light . After repeating this procedure 
on our other data-sets with /3 = 5.29, we extrapolated the results linearly in the sea quark mass to 
the physical pion mass. Combining this result with the value obtained from a similar analysis using 
6% x , we finally find 



(£)^(4GeV 2 ) = 0.0275(5)(17) , af (4GeV 2 ) = 0.0453(9)(29) . 



(3.3) 



Again, the first error is statistical and the second systematic. The magnitude of this quantity has 
been debated in the literature quite recently, see the discussion in The value we find is more 
accurate than the sum rule estimate af = 0.05(25), but fully compatible both in sign and in mag- 
nitude. The independent investigation in [J3j] finds a value of af (4GeV 2 ) = 0.053(5) which is also 
compatible with our finding. 
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Figure 3: Extrapolation of the second moment of the kaon distribution amplitude, (% 2 }k at /I 2 = 4GeV 2 , at 
/3 = 5.29, JQea = Kiigh( = 0.1359 and several values of fc va i = K he . avy . The extrapolation to the physical value is 
shown at the vertical line. The blue data point indicates the degenerate mass case, JC heavy = Ki ight . 



4. Summary 

We have performed a calculation of the first two moments of the distribution amplitudes of 
pseudoscalar mesons. In the case of the pion we have found agreement with previous sum rule 
calculations and phenomenological evaluations from experimental data. In the case of the kaon 
we have also found agreement with previous sum rule estimates, but with smaller statistical errors. 
The uncertainty due to the continuum extrapolation has been estimated and found to be small. 
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Figure 4: Interpolation of the first moment of the kaon DA, (E,)k at j8 = 5.29, rc sca = 0.1359 as obtained 
from <??4 4 . The physical point corresponding is indicated by the vertical line, the value obtained by a linear 
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